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Dynamical frictional phenomena are studied theoretically in a two-chain model with incommensurate 
structure. A perturbation theory with respect to the interchain interaction reveals the contributions 
from phonons excited in each chain to the kinetic frictional force. The validity of the theory is 
verified in the case of weak interaction by comparing with immerical simulation. The velocity and 
the interchain interaction dependences of the lattice structure are also investigated. It is shown that 
peculiar breaking of analyticity states appear, which is characteristic to the two-chain model. The 
range of the parameters in which the two-chain model is reduced to the Frenkel-Kontorova model 
is also discussed. 
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I. INTRODUCTION 

The sliding velocity dependence of the kinetic frictional 
force is of particular importance in the field of tribology. 
The Coulomb- Amonton's law states that the kinetic fric- 
tional force does not depend on the sliding velocity 0. 
It is well known that this law holds well under usual con- 
dition, but breaks in some cases. Its physical explana- 
tion is not well-established. In recent years the study of 
atomic scale frictional phenomena has been a hot issue 
|^|-^. In such microscopic systems, new physical laws 
and concepts of friction is considered to exist. Several 
works on kinetic friction have been reported using mi- 
croscopic or atomic scale models [||. Among these stud- 
ies, the Frenkel-Kontorova model (FK model) [|j and re- 
lated ones have been investigated extensively by several 
researchers . The FK model consists of the atoms 

interacting each other via harmonic force under a peri- 
odic potential. In the case that the ratio between the 
mean atomic spacing and the period of the potential is 
irrational, i.e., in an incommensurate case, the FK model 
shows a phase transition, the so-called Aubry transition 
|0,|j. When the amplitude of the periodic potential is 
smaller than a certain critical value, the maximum static 
frictional force vanishes. Above the critical amplitude 
it becomes finite. The dynamical properties of the FK 
and related models have been reported by many authors. 
Coppersmith and Fisher investigated the behavior near 
the maximum static frictional force in the FK model from 
the viewpoint of critical phenomena ||9|,|l0||. Then, the 
critical exponent of the velocity-force characteristics was 
evaluated. Sokoloff developed a perturbation theory of 
frictional force He employed a three-dimensional 

FK model, which consists of a three-dimensional de- 
formable lattice and a rigid substrate, which makes pe- 



riodic potential. He discussed frictional effects in com- 
mensurate and incommensurate systems and found that 
strong reduction of kinetic frictional force for an incom- 
mensurate system compared with that for a commen- 
surate one. The existence of a superlubricity state, in 
which the frictional force vanishes, was claimed by Hi- 
rano and Shinjo in a three-dimensional FK model Jl2| , p^ . 
Since they assumed an energy conservative system, the 
kinetic frictional force vanishes by the energy recurrence 
effect between the motion of the center of gravity and 
the lattice vibration. Matsukawa and Fukuyama inves- 
tigated the static and kinetic frictional forces based on 
a one-dimensional model ||l^,0. The model employed 
in their study consists of two atomic chains, where the 
interchain interaction, the harmonic intrachain interac- 
tion and the effect of energy dissipation are taken into 
account. In their model a definitive expression of both 
of static and kinetic frictional forces can be derived. Us- 
ing a numerical simulation, they found that the critical 
amplitude of interchain potential where the Aubry tran- 
sition occurs depends strongly on the elasticity of the 
chains. They investigated also the velocity dependence 
of the kinetic frictional force. When the interchain inter- 
action is weak, the crossover behavior of the kinetic fric- 
tional force between velocity-strengthening to velocity- 
weakening is observed as velocity increases. It was also 
found that the crossover velocity and the strength of the 
kinetic frictional force make obvious difference whether 
the lower chain is deformable or rigid. For strong inter- 
chain interaction, the velocity-strengthening behavior is 
smeared out because of large maximum static frictional 
force. Quasi-periodic and chaotic sliding states were dis- 
cussed by Elmer, Strunz and Weiss in the underdamped 
regime of the FK and FK-Tomlinson models JTot . 

For the theoretical understanding of the dynamics of 
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frictional phenomena, it is desirable to clarify the behav- 
ior and features of kinetic frictional force both by numer- 
ical and analytical methods. In this paper, using the one- 
dimensional two-chain model of friction proposed in Ref. 
iQ , we revisit the dynamic properties of the model and 
calculate the kinetic frictional force by employing both of 
perturbation theory and numerical calculation. On the 
basis of a theoretical expression of frictional force , we 
can successfully formulate the perturbation theory on the 
calculation of the kinetic frictional force for the two-chain 
model. The role of the phonon excited in each chain can 
be clarified. Numerical simulations are employed to clar- 
ify the validity of the perturbation theory and to investi- 
gate the atomic configuration and motion. The velocity 
dependence of the kinetic frictional force under various 
conditions is investigated. The physical meaning of the 
earlier results on the kinetic friction is discussed. We 
also investigate the hull function, which characterize the 
lattice structure. It is found that peculiar breaking of 
analyticity states appear, which is characteristic to the 
case of two deformable chains. Lastly we discuss the dif- 
ference between the two-chain model and the FK model 
and the parameter regime in which the two-chain model 
is reduced to the FK model approximately. 



II. TWO-CHAIN MODEL OF FRICTION 

The two-chain model of friction is reviewed in the fol- 
lowing. We consider two atomic chains, i.e., an upper 
chain and a lower chain. The atoms in each chain have 
one-dimensional degree of freedom parallel to the chain 
and interact each other via the harmonic force. Inter- 
chain interaction works between atoms in the upper chain 
and those in the lower one. The effects of energy dissi- 
pation in both chains are assumed to be proportional to 
the difference between the velocity of each atom and that 
of the center of gravity. The upper chain is subject to 
the external force parallel to the chain. Here we assume 
overdamped motion, and then the equations of motion 
of the atoms in the upper and lower chains are expressed 
as, 

rUaJaiUt - {Ui)i) = Ka{Ui+i + - 2Ui) 

+ ^Fi{u,~Vj)+Fe:,, (1) 

mbjbivi - {vi)i) = Kb{vi+i + Vi^i - 2vi) 

+ J2 Fii^'^ - "j) - Ks{v., - icb), (2) 

jea 

where Ui (vi), rua {nib), la (lb), Ka {Kb) and Na {Nb) 
are the position of the i-th atom, the atomic mass, the 
parameter of energy dissipation, the strength of the in- 
teratomic force and the number of atoms in the upper 



(lower) chain, respectively. Kg is the strength of the in- 
teratomic force between the lower chain and the substrate 
and (. . .)i represents the average with respect to i. Fj 
and Fex are the interchain force between the two atomic 
chains and the external force, respectively. We adopt the 
following interatomic potential: 
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K, 



■ exp 



Cb 



(3) 



where Ki is the strength of the interchain potential, Cb 
the mean atomic spacing of the lower chain. The inter- 
atomic force is given by Fi{x) = —-^Uj. The time- 
averaged total frictional force of the present model is 
given by the following equation. 



;ifric 



ex/t 



(4) 



It should be noted that the above expression of the fric- 
tional force is valid for both static and kinetic ones. In 
this study we apply Eq. (4) to the evaluation of kinetic 
frictional force. 



III. PERTURBATION THEORY OF FRICTIONAL 
FORCE FOR TWO-CHAIN MODEL 

We construct a perturbation theory for the kinetic fric- 
tional force that arises between two interacting atomic 
chains. In a steady state the upper chain is assumed to 
be sliding at a constant velocity, V, on average. Then 
the atomic positions in the upper and lower chains are 
expressed as. 



Ui^Vt + iCa 

Vi = icb + Svi 



6ui 



(5) 
(6) 



where Ca{cb) and 6ui{Svi) denote the mean lattice spacing 
and the deviation of atomic position from regular site 
in the upper (lower) chain, respectively. By substituting 
them into Eq. (4), we can expand the equation in terms 
of 6ui and Svi, 

Na Ni 

Y^Y.{Fi{u,-v,))t 

iGa jGb 

Na Nt 

= '^^{Fl{Vt + + SUi - jCb - SVj))t 
iGa j^b 

Na Ni 

d 



+ (^-F/(Wi - Vj)\ui=Vt+iCa,Vj=]CtSUi)t 

d 

+ {-g;^Fj{Ui - Vj)\u^=Vt+iCa,Vj=jCtSVj)t 



(7) 
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In this equation the first term vanishes because the 
present model has an incommensurate lattice structure 
and then the summations with respect to i and j cancel 
out. The second term represents the contribution from 
the phonon excited in the upper chain when the atoms 
in the lower chain are fixed at the regular sites. On the 
other hand, the third term represents the contribution 
from the phonon excited in the lower chain when the up- 
per chain is rigid, but moving at a constant velocity V 
relative to the lower chain. Then the frictional force is 
simply given by the sum of both the lowest-order terms. 
Because atomic displacements Sui and 6vi are caused by 
the interchain interaction, the calculation of the frictional 
force in the lowest order perturbation, Eq. (7), assumes 
that 6ui 7^ and Svi = in the second term, and Sui = 
and Svi ^ in the third term, respectively. The contri- 
butions from terms in which both Sui and Svi are nonzero 
arise from higher order perturbations. 

First of all, we examine the second term in Eq. (7). 
It is of particular interest because the friction in the FK 
model is given by this term as discussed just below. The 
force that acts on an atom in the upper chain in the 
leading order of Svj is expressed by the following Fourier 
series. 



Ui 



-OAlKi sin 27r 



Cb 



0.00057JOsin ( 47r--^ ) . (8) 



Hence, we can successfully neglect the second term in 
Eq. (8). Then, the equation of motion of the upper 
chain in the present perturbation theory approximately 
corresponds to the following equation of motion. 



-0.47X/sin 



2Sui} 



(Vt + lCa 

Cb 



5ui 



(9) 



This is nothing but the overdamped equation of motion 
of the FK model. 

In the perturbation theory [pT|jl^ ] , the atomic displace- 
ment is expressed by the recursive equation. 



5u, = ^^^IEl y [ dt'G.At ~ t') 
ma ^ J 



sm 



27r 

Cb 



{Vt' + iCa+5Uj) 



sm 



27r 

— {Vt + iCa + 5Ui) 
Cb 



(10) 



where Gij{t — t') is the Green function of phonons defined 
as 



k 

r2(fc)^ = (1 — cos/cCq) . 



In the lowest order, Eq. (10) is reduced to 
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(12) 
(13) 
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— iVt'+JCa) 
Cb 



(14) 



Then the kinetic frictional force per atom in the upper 
chain resulting from the phonon excitation there, F^pp°'^, 
is given in the lowest order of Kj as. 
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(15) 



As easily seen by this equation there exists the crossover 
behavior of the velocity dependence of the kinetic 
frictional force between the velocity-strengthening and 
velocity-weakening. This is the essential feature of the 
kinetic frictional force in the case of small interchain in- 
teraction in the present model, which is observed in the 
numerical simulation in Ref. [ p^ and in section V. The 
crossover velocity is determined by three characteristic 
quantities: the phonon frequency at the wavenumber cor- 
responding to the period of the underlying potential, the 
so-called washboard frequency 2TrV/cb and the damping 
constant. In order to consider the reason of this crossover 
behavior, it is interesting to see the average of the squared 
velocity fluctuation, {6ui{t)^)i^t given by. 



ma J 



( 221 



V 



(16) 



The squared velocity fluctuation shows velocity- 
strengthening behavior in the low velocity regime and 
then saturates. The former behavior comes from the 
increase of V, itself. The latter one results from the 
suppression of the fluctuation of position, {Sui{t)'^)i^t, 
in the large velocity regime, as mentioned below. We 
note here that the kinetic friction in the present model 
is caused by the energy dissipation due to the damp- 
ing of the phonons excited by the sliding motion. Hence 
the energy dissipation per unit time in the upper chain, 
^uppery^ is equal to ma-fa{SMtf)i,t- Then the fric- 
tional force, which is the energy dissipation per unit 
sfiding distance, is equal to TOa7a {Su,{t) )^,t/V. So the 
velocity-strengthening behavior of F^pp°'^ in the low veoc- 
ity regime comes from that of {5ui{t)^)i,t- As V becomes 
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greater than (f )^ (Sudtfk. saturates and then 

the frictional force decreases. 

Next we consider the kinetic friction resulting from the 
phonon excitation in the lower chain. In this case we can 
express the interatomic force that act on an atom in the 
lower chain using a Fourier series as follows. 



- Vt 

-0.83X/sin 2ii^ 



-0.098K/sin At: 



Vi - Vt 



(17) 



We can neglect the second term. The strength of inter- 
chain interaction is not equal to that in Eq. (8). This 
comes from the spatial profile of the interchain potential 
and the summation on the upper lattice site. Then the 
equation of motion is given by 



"27r 



2n 

'OMKi — cos 



-i-Vt + iCb) 



Sv,. (18) 



After a similar calculation to that of F"pp", we get the 
kinetic frictional force resulting from the phonon excita- 
tion in the lower chain. It is written as ^1°^°'' and is 
given by 



;ilowcr 



Nb {OMKj 



V 



Na 2^bmb 



(19) 



where 51b (fc) — \/ {Kg + 2Kb (1 — cos kcb)) /mb is the 
phonon frequency of the lower lattice. ^1°™'' is also eval- 
uated per atom in the upper chain. 

The total kinetic frictional force F^'^^ is given as — 
^uppor flower r^^Ycn the Strength and velocity depen- 
dence of the total kinetic frictional force are determined 
by both of F"PP" and i^i°wcr xhis is the essential differ- 
ence with the FK model. In the latter F^""- = i^"PP<='-, of 
course. 

The perturbation theory is expected to be valid under 
the condition that the atomic displacement of the lattices 
is small compared to the mean lattice spacing. That 
condition is expressed as follows. 



{5u,{tf 



0A7Ki 
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These inequalities depend on velocity because atomic dis- 
placements are suppressed by large velocity. The rea- 
son is that in the equation of motion, Eqs. (1) and (2), 
the effect of the high-velocity-sliding motion of atoms 
effectively weakens that of interchain interaction, which 
causes the atomic displacement. It is a dynamical effect 
in sliding. Here we make some comments on these condi- 
tions required for the two chains. For the perturbation in 
the FK model, only one inequality, Eq. (20), is required, 
but in the two-chain model both inequalities must be sat- 
isfied simultaneously. It is to be noted that in the present 
lowest order calculation the total kinetic frictional force 
is the sum of the contributions of phonons in the upper 
and lower chains as mentioned above. If the conditions 
Eqs. (20) and (21) are broken, there appears the cross 
term of Sui and 5vi, which expresses the contribution 
from the interaction of phonons of both chains. 



IV. METHOD OF NUMERICAL SIMULATION 

We perform a numerical calculation of the kinetic fric- 
tional force and examine the validity of the perturbation 
theory developed in the previous section. The equations 
of motion are numerically solved using the Runge-Kutta 
formula at the fourth order. We assume periodic bound- 
ary conditions in both chains. Hence the ratio of the 
mean lattice spacings of the upper and lower chains Ca/cb 
is equal to the ratio Nt/Na- We choose the ratio by using 
the continued-fraction expansion of the golden mean. 



and 



V2ma 



OMKi 



■€.Ca, (20) 



V2', 



mb 



Ca 
Cb 



Nb_ 

Na 



144 
"89" 



1.617- 



(22) 



where the expansion is truncated at the tenth order. We 
set the values of the model parameters as 



Na 
Ka 



89, Nb = 144, Ca 



144 
"89" 



, Cb = l,TOa = mfc = 1, 



I, Kb = 0,7a = 76 1, 



(23) 



where the intrachain force is neglected, i.e., Kb = 0, for 
simplicity. Such simplification was employed also in the 
earlier study ||lj]. The lower chain, therefore, consists 
of individually oscillating atoms, that is, the Einstein's 
model of lattice vibration is adopted in the lower lattice. 

Throughout this study, the frictional force is evaluated 
per atom in the upper chain. 



<Cb. (21) 



V. COMPARISON BETWEEN THE RESULTS OF 
THE PERTURBATION THEORY AND 
NUMERICAL CALCULATION 

In Section III we have presented the kinetic frictional 
force calculated in the lowest order perturbation. Higher 
order terms of the perturbation can be calculated in a 
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similar manner, but it is easier and more straightforward 
here to check the vahdity of the lowest order perturba- 
tion by comparing its result with that of numerical cal- 
culation. First of all, we must investigate the relation- 
ship between the strength of the interchain interaction 
and the validity of the perturbation theory in the case 
that the lower chain is fixed, which corresponds to the 
FK model as mentioned in Section III. In this case the 
contribution to the kinetic frictional force results only 
from the phonon excitation in the upper chain. Even in 
such a single deformable chain case the validity of the 
perturbation theory has never been examined so far es- 
pecially for strong interchain interaction. The kinetic 
frictional force in the case of weak interchain interac- 
tion, Kj ~ 0.1, which is less than the critical value of 
the Aubry transition, is shown in Fig. 1(a) as a func- 
tion of the sliding velocity. Then the maximum static 
frictional force vanishes because of the absence of the 
Aubry transition. A good agreement between the re- 
sult of the lowest order perturbation theory and that of 
the numerical calculation is obtained and the velocity- 
strengthening and velocity- weakening features are clearly 
observed. From this result, it is confirmed that the re- 
sult of the lowest order perturbation theory is quite good. 
This is consistent with the condition, Eq. (20), which 

gives \J {5ui{t)^)i^t ~ 0.01 ^ Ca even in the case of van- 
ishing velocity. Hence it is found that dominant con- 
tribution to the kinetic frictional force comes from the 
lowest order, and then higher order terms are negligi- 
ble. On the other hand, for strong interchain interac- 
tion, Kj — 1, which is greater than the critical value of 
the Aubry transition and therefore finite static frictional 
force exists as observed in Fig. 1(b), obvious discrepancy 
is observed between the result of the perturbation theory 
and that of the numerical calculation in the low velocity 
regime. This dicrepancy appears to arise from the Aubry 
transition. The effect of the Aubry transition affects seri- 
ously the kinetic friction in the low velocity sliding state. 
However, we should note that in the high velocity regime 
{V > 0.6) the result of the lowest order perturbation and 
that of the numerical calculation agree well. This is also 
consistent with Eq. (20). The solid lines in Fig. 1 denote 

{Sui{t)'^)i^t/ca- As seen from Fig. 1(b), the frictional 
force derived from the lowest order perturbation theory 
agrees well with that from numerical calculation, where 

the velocity is so large that ^J {5ui{t)^)i^t ^ Ca is satis- 
fied. Therefore, in this velocity regime, it is considered 
that the sliding state is well described by the lowest or- 
der perturbation theory. In the following, we analyze the 
lattice structure of the sliding chain and clarify that the 
discrete features of the atomic distribution caused by the 
Aubry transition is relevant to the above discrepancy in 
the low velocity regime. 

In order to show that the validity of the perturbation 



theory has a close relationship to the lattice structure of 
the sliding chain, it would be useful to consider the hull 
function, h{x). It is defined as follows, 

Ui = iCa + h{ica). (24) 

The hull function is a periodic function, the period of 
which is that of the potential: h{x) — h{x -\- Cb), and 
therefore the variable of the hull function is defined in the 
range from to C{,. We pay our attention first to the static 
case {V — 0). When the strength of the interchain inter- 
action is less than the critical value of the Aubry tran- 
sition, the hull function is continuous and almost sinu- 
soidal. This reflects the sinusoidal form of the interchain 
force of Eq. (8). Otherwise, however, the hull function is 
discontinuous and has gaps at several positions because 
of the breaking of analyticity due to the Aubry transi- 
tion as shown in Fig. 2(a). In particular, among the 
gaps, the largest gap is located at the half of the period, 
Cf,/2. This largest central gap is considered to charac- 
terize the Aubry transition and means that the atomic 
distribution vanishes at the potential maximums and the 
atoms are bound strongly near the potential minimums. 
We show here that the hull function gives us an interest- 
ing insight about the atomic motion in sliding under the 
influence of the underlying potential as noticed in Ref. 
Jigf . Using such a dynamical hull function, the atomic 
position in sliding is given by Ui — ICa + Vt + h{ica + Vt). 
Figs. 2(b)~(e) show some snapshots of the hull function 
at various steady sliding velocities for the strong inter- 
chain interaction, Kj = 1. When the velocity is small 
enough (Fig. 2(b)), the reminder of the gaps is clearly 
observed and the shape is highly distorted. It is due to 
the large lattice distortion caused by the Aubry transi- 
tion, which is responsible for finite static frictional force. 
The lattice still remains highly distorted even in a slid- 
ing state. Such sliding state is characterized by the al- 
most discontinuous and nonuniform spatial distribution 
of atoms which is caused by local atomic stick-slip motion 
in which most atoms stay nearby the potential minimums 
and other atoms jump over the potential maximums. It 
should be noted that such slip motion, of course, needs 
finite time and then the hull function has no gap in the 
case of finite velocity. In Fig. 2(b) the reminders of some 
small gaps are vanishing, but that of the largest gap is 
apparent. The deviation of the position of the reminder 
of the largest gap from the half of the period is due to 
the external force added to the sinusoidal interchain force 
pof . Because the effect of large lattice distortion due to 
the Aubry transition is not taken into consideration in 
the present perturbation theory, it fails to explain the be- 
havior of the total kinetic frictional force F^"^'^ obtained 
by the simulation. As the velocity increases, the ampli- 
tude of the hull function decreases and its form becomes 
smoother and comes close to the sinusoidal one. This 
means that the large lattice distortion is suppressed in 
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the high velocity regime. In this regime the perturba- 
tion theory gives F^"'^ which agrees well with that by the 
simulation. As mentioned before, the suppression effect 
can be confirmed directly from Eq. (20), as shown by the 
solid lines in Fig. 1. Even when the lattice structure is 
subject to the Aubry transition, the dynamic suppression 
effect washes out it in the high velocity regime and leads 
to the sliding state that is described by the lowest order 
perturbation theory. 

It turns out that the frictional property of the single 
deformable chain on a periodic potential can be under- 
stood on the basis of the lowest order perturbation the- 
ory. We move on to the two-chain model. We clarify 
the validity of the lowest order perturbation theory also 
here and where is the difference between the two-chain 
and the FK models. Figs. 3(a) and (b) show the kinetic 
frictional force for Kg = 2 and 10 respectively in the case 
of weak interchain interaction, Kj = 0.1. In these cases 
the maximum static frictional force is vanishing and the 
velocity-strengthening and velocity- weakening features of 
the kinetic frictional force are observed. The contribu- 
tions to the kinetic frictional force from the phonon exci- 
tations in the upper and lower chains, F"pp*'' and F^°'^^'^, 
and the total kinetic frictional force, F^''^'^, calculated by 
the perturbation theory are also plotted in these figures. 
For Ks = 2 (Fig. 3(a)), the contribution to the kinetic 
frictional force from the lower chain is dominant in the 
overall velocity regime. On the other hand, when the 
lower chain becomes stiffer. Kg = 10 (Fig. 3(b)), the con- 
tributions from both chains become important. Because 
the crossover velocity between velocity-strengthening and 
velocity-weakening features is different for each chain in 
this case, the velocity dependence of the total kinetic 
frictional force is modified in a velocity regime between 
the two crossover points (0.7 < V < 4). In this veloc- 
ity regime, the velocity-weakening behavior of the kinetic 
frictional force due to the lower chain and the velocity- 
strengthening one due to the upper chain arc complemcn- 
tal each other, and therefore the total kinetic frictional 
force becomes velocity-insensitive. It is found in both 
the cases that the results of the lowest order perturba- 
tion theory can explain excellently the behavior of the 
kinetic frictional force observed in the numerical simu- 
lations when the strength of the interchain interaction 
is less than the critical value of the Aubry transition. 
Therefore, it is confirmed also in the two-chain model 
that the lowest order perturbation is valid definitively in 
the absence of the Aubry transition, and then higher or- 
der contributions are quite negligible. Thus the lowest 
order perturbation gives us a correct and quite simple 
viewpoint on the kinetic frictional force of the two-chain 
model in the absence of the Aubry transition that the be- 
havior of each chain can be understood as an individual 
deformable chain on the periodic potential and then the 
kinetic frictional force can be also composed of two con- 
tributions that come from such individual two chains. If 



the Aubry transition exists in the two-chain model, such 
a simple viewpoint on two chains is modified largely and 
altered to a complicated one. Fig. 3(c) shows the kinetic 
frictional force for Kg = 2 in the case of strong inter- 
chain interaction, Ki = 1, which is the same strength 
as that in Fig. 1(b). There exists a obvious discrepancy 
between the result of the perturbation and that of the nu- 
merical calculation, but in a high velocity regime a good 
agreement is obtained. Such behavior is quite similar to 
that for the FK model in Fig. 1(b). Then we investigate 
the difference in the kinetic frictional force and the lat- 
tice structure between the FK model and the two-chain 
model. In particular, it is quite important to understand 
how the strong interchain interaction leads to the Aubry 
transition in two chains, and how both lattices behave in 
sliding states. In order to discuss lattice structures of two 
chains, we define two hull functions for the two chains of 
the present model in the following equations. 

Ui = iCa + Vt + ha{iCa + Vt), (25) 

Vi = ich -\- hb{icb). (26) 

In this case the periodicity of the hull functions is ex- 
pressed as 

ha{x) = ha{x + Cb), hb{x) = hb{x + Ca). (27) 

The sliding states of the strongly interacting two-chain 
model can also be analyzed using these two hull func- 
tions. Figs. 4(a)~(e) show the hull functions given by 
snapshots taken in the case of no external force and in 
steady sliding states at several velocities for the system 
with Kj = 1. For these model parameters, it should be 
noted that the Aubry transition occurs in the lower chain 
in the stationary state in the sense mentioned just below. 
It is confirmed from the fact that the hull function hb in 
Fig. 4(a) shows the largest gap at the half of the period, 
Ca/2 ss 0.81. On the other hand, the hull fimction for 
the upper chain is also much affected by the strong in- 
terchain interaction and have gap structure. But no gap 
exists at the half of the period, Cb/2 = 0.5. The state of 
the upper chain would be a sort of breaking of analyticity 
state, but we can distinguish that from the conventional 
Aubry's breaking of analyticity state for the FK model, 
which is characterized by the largest gap at the half of 
the period. The interchain force which affects the upper 
chain is not simple and spatially modulated because the 
lower chain is softer in elasticity than the upper one and 
is highly distorted. The atoms of the soft lower chain are 
bound near the minimum points of the almost periodic 
potential caused by the stiffer upper chain. This leads to 
the conventional Aubry's breaking of analyticity state in 
the lower chain. Thus, in the stationary state, the elastic- 
ities of two chains cause a peculiar breaking of analyticity 
state in each chain when the strength of the interchain 
interaction is greater than a critical value. The critical 
strength of the interchain interaction itself varies with 
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the elastic parameters as observed in Ref. These 
features of the hull functions are quite complicated and 
the perturbational treatment is apparently invalid in this 
case. 

In sliding states with very low velocity, both the hull 
functions, ha and hb, retain the reminder of the gap struc- 
ture and are highly distorted (Fig. 4(b)). As velocity 
increases, however, they gradually change their forms, 
approach to sinusoidal ones and reduce the amplitudes 
(Fig. 4(c)). Such behavior of the hull functions is es- 
sentially the same as that for the FK model. The sinu- 
soidal form of the hull function /if, is almost recovered at 
V ^ 0.658 as seen in Fig. 4(d), but the hull function 
ha is still distorted from the sinusoidal form. This dif- 
ference in behavior between the two hull functions would 
come from the difference of their elasticities. For a higher 
velocity (Fig. 4(e)), both the hull functions become si- 
nusoidal obviously. It is to be noted here that the main 
contribution to the kinetic frictional force comes from the 
lower chain, in which the conventional Aubry transition 
takes place as observed in Fig. 3(c). Hence the result of 
the numerical calculation on the kinetic frictional force 
comes closer to that of the perturbation theory above the 
velocity corresponding to Fig. 4(d), where h^ recovers si- 
nusoidal form. Above the velocity corresponding to Fig. 
4(e), where both hull functions become sinusoidal, they 
agree quite well. This means that the result of the nu- 
merical calculation for the two-chain model agrees well 
with that of the perturbation theory only when the ef- 
fects due to the breaking of analyticity vanish in both 
chains and therefore the sinusoidal forms appear in both 
the hull functions. 

Now we consider the conditions required for the valid- 
ity of the lowest perturbation theory, Eqs. (20) and (21). 
Fig. 5 shows ^{{Su,)^),/ca (Eq. (20)) and ^/ {(Sv,)^) , / Cb 
(Eq. (21)) as a function of velocity, where the parame- 
ters are the same as those in Fig. 4. In the small ve- 
locity regime {{5uiY)i/ Ca is quite small (« 0.06), but 
yj {{5viY) i / Cb is rather large [k, 0.3). The large atomic 
displacement of the lower chain that is obtained by the 
lowest perturbation theory means that the theory itself is 
not valid although the atomic displacement of the upper 
chain calculated with the theory is quite small. In the 
high velocity regime above the characteristic velocities 
for two chains, which is about unity, the atomic displace- 
ments of both chains are suppressed and then the lowest 
order perturbation theory gives well description of the 
present system. 

We discuss here the range of the parameters in which 
the present two-chain model is reduced to the FK model, 
i.e., a model that consists of a one-chain on a rigid peri- 
odic potential. This offers an interesting insight into the 
two-chain model because the importance of being the two 
chains is understood from it. In order to investigate the 
problem, we note the change of the gap structures of the 



hull functions. Fig. 6 shows the two hull functions, ha 
and hb, in the absence of the external force for several 
Ks's [Ki = 1 and Kb = 0). For smaU Ks {K^ = 3, Fig. 
6 (a)), as observed also in Fig. 4(a), the central gap exists 
only in hb, i.e., the conventional Aubry transition occurs 
in the lower chain. In this case the elasticity of the lower 
chain has great importance obviously because the lower 
chain has a dominat contribution to the kinetic frictional 
force. As Kg is increased {Kg — 5, Fig. 6 (b)), however, 
the central gap disappears in hb- When Kg is increased 
further {Kg > 10, Figs. 6 (c) and (d)), the central gap 
appears in ha- This is obviously contrary to those for 
small Kg^s, but this gap structure of ha is just the same 
as that shown in Fig. 2 (a) for the Aubry transition in 
the FK model. Then the amplitude of hb is strongly sup- 
pressed and all the gaps shrink. In this case the lower 
chain is very stiff, and its contribution to the frictional 
force is quite smaller than that from the upper chain in 
the low velocity regime. In the high velocity regime it is 
possible that the high- velocity-sliding upper chain excites 
high-frequency phonons in the stiff lower chain, and then 
a finite contribution to the kinetic frictional force from 
the lower chain is observable if the contribution from the 
upper chain is sufficiently small in this velocity regime. 
This is, of course, an essentially different point between 
the FK and two-chain models. If our discussion is re- 
stricted to the low velocity regime, the two-chain model 
can be reduced to the FK model in the large Kg regime. 
Kg ^10 for Kj — 1. Then, the maximum static and 
kinetic frictional forces become close to that for the FK 
model. Such behavior of the maximum static frictional 
force has been reported in Ref. jl4j. The kinetic fric- 
tional force in such a case with the same values of the 
parameters with Fig. 3 (c), Ki — 1 and Kb = 0, but Kg, 
which is equal to 20, is shown in Fig. 7, where the lower 
chain is very stiff. The magnitude of the kinetic frictional 
force in the low velocity regime comes closer to that for 
the FK model in comparison with Fig. 3 (c). In the high 
velocity regime, however, the contribution from the stiff 
lower chain nevertheless becomes important. In Fig. 8 
we summarize the characteristic parameter regimes for 
the two-chain model using the parameters, Kg and Kj. 
In the regime I, ha shows the largest central gap, while 
there is no central gap in hb, as observed in Figs. 6 (c) 
and (d). That is, as discussed above, the two-chain model 
is approximated by the FK model especially in the sta- 
tionary and the low-velocity sliding states. In the regime 
II, however, no central gap exists in ha and the gap struc- 
tures of both the two hull functions become complicated. 
Then the contribution from the lower chain to the fric- 
tional force becomes important. The values of Kg at the 
boundary between the regimes I and II increase linearly 
with Kj. It is to be noted that there would be no exact 
phase boundary between the regimes I and II. 

Thus the elasticity of two chains and the interchain in- 
teraction for the two-chain model affect each other com- 
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plexly when the strength of the interchain interaction 
is so strong that the Aubry transition occurs. Then 
the lowest order perturbation theory is not appUcable to 
two-chain models in that case. Furthermore the feature 
characteristic of the two-chain model is more complicated 
than that for the FK model. It is found, however, that 
for both models the validity of the lowest order perturba- 
tion theory is retrieved in the high velocity regime, and 
the two-chain model is reduced to the FK model when 
the lower chain becomes stiffcr than the upper chain un- 
der a certain condition on the strength of the interchain 
interaction. 

VI. SUMMARY 

We have investigated the kinetic frictional force for a 
one-dimensional two-chain model with an incommensu- 
rate lattice structure. On the basis of the explicit theo- 
retical expression of frictional force, we have formulated 
a perturbation theory of the kinetic frictional force and 
clarified the contribution from the phonon excitation in 
each chain to the kinetic frictional force. From the study 
on the FK model (single deformable chain model) and 
the two-chain model we have found that the lowest or- 
der perturbation theory explains the numerical results 
excellently over a wide range of velocity if the strength 
of the interchain interaction is less than the critical value 
of the Aubry transition. Even if the finite static frictional 
force due to the Aubry transition appears, the lowest or- 
der perturbation theory is still valid in the high velocity 
regime where the suppression of atomic displacement be- 
comes significant and the large lattice distortion due to 
the breaking of analyticity is reduced by high-velocity 
sliding motion. In particular, for the two-chain model 
with certain elastic parameters, we have found that the 
conventional Aubry's breaking of analyticity state, which 
is observed in the FK model and is characterized by the 
largest central gap of the hull function, and an another 
breaking of analyticity state appear in the lower and up- 
per chains, respectively. The latter lacks the largest cen- 
tral gap and is not well-defined in the context of the 
conventional Aubry transition. These stationary states 
cause complicated sliding lattice structures which are re- 
sponsible for the discrepancy between the result of the 
perturbation theory and that of the numerical simula- 
tion. As velocity is increased, however, these states of 
two chains come to sliding states described by the lowest 
order perturbation theory. Thus, because each chain play 
an important role on the kinetic friction, it is necessary 
to take into account the contributions to the kinetic fric- 
tional force from both chains for understanding kinetic 
frictional force of the two-chain model |2l] . 

We have also clarified the range of the parameters in 
which the two-chain model is reduced to the FK model. 
In the suitable range of Kj and Kg the former reduces to 



the latter in the low-velocity regime. Even in that range, 
however, discrepancy between them appears in the high- 
velocity regime. 

The present results for the two-chain model is use- 
ful to understand the earlier results in Refs. In 
these studies it was found numerically that the kinetic 
frictional force shows crossover behavior from velocity- 
strengthening to velocity-weakening and the crossover 
velocity and the strength of the kinetic frictional force 
are obviously different between the FK and two-chain 
models. Their results can be understood well with the 
perturbation theory of the present study. 

In the present study the parameters of the model are 
chosen artificially and are not necessarily considered the 
correspondence to those in realistic materials because the 
system is one-dimensional. Actual frictional phenomena, 
however, take place in two- or three-dimensional systems. 
The effect of the three-dimensionality of phonon excita- 
tion and the elasticity in each substance would be cru- 
cial to the kinetic friction in such systems. Moreover, 
although the present model has no randomness, random- 
ness such as impurities would affect the kinetic and static 
frictional forces. The relationship between the velocity 
dependence of the kinetic frictional force and the static 
frictional force caused by pinning due to such random- 
ness is not well established. These issues will be discussed 
elsewhere p^ . 
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Fig. 1 

Kinetic frictional force and averaged atomic displacement 
plotted against velocity (for the FK model), (a) A weak 
interchain interaction case: Ki = 0.1. (b) A strong in- 
terchain interaction case: Kj = 1. Marked and dotted 
lines represent the result of the numerical simulation and 
that of the perturbation theory, respectively. Solid line 
without marks denotes the averaged atomic displacement 
calculated with Eq. (20), which is renormalized by the 
mean lattice spacing. 

Fig. 2 

Hull functions for the FK model with strong interchain 
interaction Kj = 1. The velocities correspond to the fig- 
ures are (a) (no external force), (b) 2.03 x 10~^, (c) 
1.11 X 10-\ (d) 3.95 X 10-1, and (e) 3.22. 



Fig. 3 

Kinetic frictional force plotted against velocity (for the 
two-chain model), (a) and (b) correspond to a weak in- 
terchain interaction case Kj = 0.1, where only Kg is 
different: (a) Kg = 2 and (b) Kg = 10. (c) is for strong 
interchain interaction Kj = 1 and Kg = 2. Marked line 
represents the results of the numerical simulation. Solid 
line is the kinetic frictional force obtained by the per- 
turbation theory, which is decomposed into the contri- 
bution from the upper chain (dotted line) and that from 
the lower chain (broken line). 

Fig. 4 

Hull functions for the two-chain model with strong in- 
terchain interaction {Kj = 1) and Kg = 2. The graphs 
in the left (right) row are the hull functions for the up- 
per(lower) chain. The velocities are (a) (no external 
force), (b) 1.36 x 10"^ (c) 1.52 x lO'S (d) 6.58 x lO'S 
and (e) 1.33. 

Fig. 5 

Averaged atomic displacements plotted against velocity. 
The parameters are the same as those in Fig. 4. Solid 
and dashed lines represent ^/ {{5ui)'^)i/ Ca calculated with 
Eq. (20) and \/ {{5viY)i/ Cb calculated with Eq. (21), re- 
spectively. 

Fig. 6 

Hull functions for the two-chain model with strong inter- 
chain interaction Ki = 1. /^g's chosen are (a) 3, (b) 5, 
(c) 10, and (d) 16. 

Fig. 7 

Kinetic frictional force plotted against velocity. Marked 
line with closed squares represents the results of the nu- 
merical simulation for Kj = 1 and Kg = 20. Marked line 
with open circles represents the results of the numerical 
simulation for the FK model with Ki = 1. Solid line is 
the kinetic frictional force obtained by the perturbation 
theory, which is decomposed into the contribution from 
the upper chain (dotted line) and that from the lower 
chain (broken line). 

Fig. 8 

Characteristic parameter regimes for the two-chain 
model. In the regime I, ha shows the largest central gap, 
which is essentially the same as that for the FK model, 
and all the gaps of hi, highly shrink and no central gap 
exists there. In the regime II, no central gap exists in ha, 
and both ha and hi, show complicated gap structures. 
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